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Abstract 

The  heat  and  mass  transfer  characteristics  in  a  steam  reformer  are  investigated  via  experimental  and  numerical  approaches  and  a  new  configuration 
of  packed  catalysts  is  proposed  for  effective  hydrogen  production.  Prior  to  the  numerical  investigation,  parametric  studies  are  carried  for  the  furnace 
temperature,  steam-to-carbon  (S:C)  ratio,  and  gas  flow  rate.  After  validation  of  the  developed  code,  numerical  work  is  undertaken  to  determine  the 
relationship  of  the  operating  parameters.  Based  on  the  experimental  and  numerical  results,  and  with  the  goal  of  obtaining  optimum  heat  transfer 
characteristics  and  an  efficient  catalyst  array,  a  new  configuration  for  the  packed  bed  is  proposed  and  numerically  investigated  taking  into  account 
the  endothermicity  of  the  steam  reforming  reaction.  A  bed  packed  repeatedly  with  inert  and  active  catalysts  is  found  to  be  an  efficient  means  to 
obtain  the  same,  or  better,  hydrogen  production  with  small  amounts  of  the  active  catalysts  compared  with  a  typical  steam  reformer. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Offering  high  energy-conversion  efficiency  and  low  emis¬ 
sion  of  air  pollutants,  the  fuel  cell  holds  promise  for  application 
as  an  energy  source.  The  various  types  of  fuel  cell  developed 
to  date  are  operated  with  hydrogen  [1-3].  For  successful  com¬ 
mercialization,  a  stable  supply  of  hydrogen  is  required  at  low 
cost  and  with  high  efficiency  [4].  There  are  many  technolo¬ 
gies  to  produce  hydrogen  from  hydrocarbons  [5-7].  Among 
these  technologies,  industrial  steam  reforming  (SR)  reactors 
with  various  types  of  burner  have  been  widely  studied.  Although 
they  have  a  strong  endothermic  nature,  SR  reactors  deliver  a 
high  hydrogen  yield  [4,8].  There  are  also  various  technologies 
to  supply  the  necessary  reaction  heat  in  an  efficient  man¬ 
ner. 

Basic  SR  reactor  designs  vary  according  to  the  arrangement 
of  the  burners,  namely  top-fired,  bottom-fired,  side-fired  and  ter¬ 
race  wall  types.  Different  tube  wall  temperatures  and  heat  flux 
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profiles  can  be  observed  according  to  the  position  of  the  burn¬ 
ers  [9].  Tubes  packed  with  supported  nickel  catalysts  should  be 
selected  rigorously  due  to  material  issues  related  to  exposure  at 
high  temperatures.  In  addition,  the  large  amount  of  heat  sup¬ 
ply  and  the  endothermic  reforming  reactions  cause  significant 
axial  and  radial  temperature  gradients  in  the  reformer.  These 
phenomena  affect  the  performance  of  the  reformer,  especially 
in  the  central  region,  because  of  heat  transfer  limitations.  Since 
combustion  is  controlled  via  dilution  with  excess  air,  the  mate¬ 
rial  problem  can  be  solved.  Basically,  the  temperature  at  the 
tube  wall  should  be  maintained  at  a  value  that  is  as  low  as 
possible  in  order  to  extend  the  tube  lifetime  [10,11].  To  com¬ 
pensate  for  shortcomings  such  as  reduction  in  the  lifetime  of  the 
tube,  various  approaches  have  been  reported,  including  distri¬ 
bution  of  active  catalysts  and  optimization  of  heat  flux  profiles 
[12,13]. 

Reactor  analysis  and  design  are  necessary  in  order  to 
develop  different  catalyst  arrangements  for  high  efficiency 
and  long-term  stability  with  optimized  heat  transfer.  To  this 
end,  a  Langmuir-Hinshelwood  type  of  heterogeneous  reac¬ 
tion  model  may  be  adequate  for  analyzing  the  heat  and 
mass  transfer  characteristics  between  bulk  gas  and  catalysts, 
instead  of  the  pseudo-homogeneous  model  [14].  In  addition, 
an  effectiveness  factor  could  be  used  to  take  into  account 


S.  Lee  et  al.  /  Journal  of  Power  Sources  180  (2008)  506-515 


507 


Nomenclature 

as f  interfacial  surface  area  (m-1) 

Cp  heat  capacity  ( J  kg- 1  K- 1 ) 

D  diffusivity  (m2  s  - 1 ) 

E  activation  energy  (kJ  kmol- 1 ) 

h  heat  transfer  coefficient/mass  transfer  coefficient 

(W m-2  K-1)  or  (ms-1) 
k  thermal  conductivity  ( W  m- 1  s  - 1 ) 

K  permeability  (m2) 

p  pressure  (N  m-2) 

q"  heat  flux  (W  m-2) 

r  radial  coordinate  of  the  reactor  (m) 

R  Universal  gas  constant  (8 . 3 1 4  J  mol- 1  K- 1 ) 

t  time  (s) 

T  temperature  (K) 

u  axial  velocity  (ms-1) 

v  radial  velocity  (ms-1) 

W  molecular  weight  (kg  kmol - 1 ) 

Y  mass  fraction 

z  axial  coordinate  of  the  reactor  (m) 

Greek  symbols 
s  porosity 

rj  effectiveness  factor 

pi  dynamic  viscosity  (kg  s- 1  m- 1 ) 

p  density  (kg  m- 3 ) 

Superscript 
in  inlet 

Subscripts 

c  centerline  of  the  reactor 

cat  catalyst 

d  mass 

D  Darcian 

eff  effective 

f  fluid  phase 

s  solid  phase 

t  heat 

w  wall 


2.  Experiments 

2.7.  Experimental  set-up 

A  schematic  diagram  of  the  SR  experimental  set-up  is  pre¬ 
sented  in  Fig.  1.  A  fixed-bed  reactor  is  located  in  an  electrical 
furnace.  Air,  fuel  and  water  are  fed  as  a  mixture  into  the  reac¬ 
tor  while  the  temperature  is  controlled  by  a  furnace.  Water  is 
vapourized  by  an  external  heat-exchanger  and  supplied  to  the 
reactor.  The  reactor  is  a  stainless-steel  tubular  bed  (inner  diame¬ 
ter  =  15.7  mm,  length  =126  mm)  that  is  filled  with  nickel  catalyst 
supported  on  alumina;  the  tube  is  0.5-mm  thick.  Reactor  temper¬ 
atures  are  measured  and  monitored  continuously  at  eight  points 
along  its  centre  line  by  means  of  K- type  thermocouples.  Temper¬ 
atures  at  the  outer  wall  of  the  reactor  are  also  measured  at  several 
locations.  Pipelines  that  connect  the  evaporator  and  the  reactor 
are  heated  electrically  to  prevent  steam  condensation.  The  flow 
rates  of  methane,  nitrogen  and  air  are  controlled  by  mass-flow 
controllers  (MFCs),  and  the  flow  rate  of  water  is  controlled  by  a 
high-performance  liquid  chromatograph  (HPLC)  pump.  Mois¬ 
ture  in  the  product  gases  is  removed  by  a  chiller  and  a  silica  gel 
trap  for  analyses  of  gaseous  compositions.  Dry  product  gases  are 
analyzed  with  an  Agilent  6890  gas  chromatograph  (GC)  using 
a  thermal  conductivity  detector  (TCD)  and  a  flame  ionization 
detector  (FID).  FU,  CO,  CO2,  and  N2  are  detected  by  the  TCD 
with  argon  as  a  carrier  gas,  and  all  hydrocarbons  are  detected  by 
the  FID  with  helium  as  a  carrier  gas. 

2.2.  Experimental  results 

In  the  present  experiment,  the  influence  of  the  furnace  tem¬ 
perature  and  the  steam- to-carbon  (S:C)  ratio  are  assessed  with 
an  average  catalyst  pellet  size  of  250-425  jim  to  exclude  the 
effect  of  catalyst  shape.  Additionally,  heat  transfer  and  diffusion 
limitations  are  studied  with  different  volume  flow  rates  of  the 
feedstock.  To  take  the  diffusion  limitation  into  account,  3 -mm 
spherical  catalysts  are  used. 

In  order  to  investigate  the  characteristics  of  commercial  cata¬ 
lysts,  the  furnace  temperature  is  varied.  Fig.  2  shows  the  variation 
in  the  composition  of  the  production  gas  with  furnace  tempera- 


the  significant  reduction  of  reaction  kinetics  due  to  intra¬ 
particle  mass-transport  limitations  [15].  In  this  study,  heat  and 
mass  transfer  phenomena  are  extensively  analyzed  using  a 
heterogeneous  two-dimensional  model  with  experimental  and 
numerical  approaches.  Due  to  heat  transfer  characteristics  and 
low  residence  time  for  chemical  species  participating  in  the 
reactions,  hydrogen  production  at  high  reactant  flow  rate  is 
remarkably  decreased.  Improved  configurations  of  the  catalytic 
bed  that  entail  repeated  inert  and  active  catalyst  packing  are 
introduced  with  the  aim  of  mitigating  the  heat  transfer  limita¬ 
tion. 


Fig.  1.  Schematic  of  experimental  setup. 
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Fig.  2.  Product  gas  concentrations  vs.  furnace  temperature  (S:C  =  3, 
GHSV  =  5000h-1). 


ture.  The  gas  hourly  space  velocity  (GHSV)  is  defined  as 


GHSV  (h-1)  = 


volume  flow  rate  of  feedstock 
volume  of  catalyst  bed 


(1) 


In  this  case,  GHSV  is  5000 h-1  and  the  S:C  ratio  is  fixed 
at  3.0.  In  Fig.  2,  hydrogen  production  tends  to  be  higher  and 
follows  the  equilibrium  composition  of  the  species  as  the  tem¬ 
perature  is  increased.  High  temperature  is  favourable  for  the 
reaction,  but  energy  input  is  proportional  to  the  reaction  tem¬ 
perature.  Based  on  the  thermodynamic  equilibrium,  hydrogen 
production  is  highest  at  around  700  °C  under  the  given  operat¬ 
ing  conditions.  The  variance  in  hydrogen  production,  however, 
is  negligible  when  the  temperature  is  above  700  °C. 

The  influence  of  the  S:C  ratio  was  investigated  at  850  °C,  at 
which  GHSV  is  5000  h_1 .  As  expected,  a  higher  S:C  ratio  results 
in  higher  hydrogen  yield,  as  shown  in  Fig.  3.  The  temperature  is 
sufficient  to  give  a  higher  hydrogen  yield  with  respect  to  injected 
methane. 


Fig.  3.  Product  gas  concentrations  vs.  S:C  ratio  (furnace  temperature  =  850  °C, 
GHSV  =  5000h_1). 


Fig.  4.  (a)  Temperature  distribution  of  catalyst  and  outer  wall  (furnace  tem¬ 
perature  =700°  C,  GHSV  =  2500 h-1,  S:C  =  3).  (b)  Temperature  distribution 
of  catalyst  and  outer  wall  (furnace  temperature  =  700  °C,  GHSV  =  5000  h-1, 
S:C  =  3).  (c)  Temperature  distribution  of  catalyst  and  outer  wall  (furnace  tem¬ 
perature  =  750  °C,  GHSV  =  10,000 h”1 ,  S:C  =  3). 


Finally,  in  order  to  determine  the  effect  of  heat  trans¬ 
fer,  a  spherical  nickel  alumina  catalyst  with  a  size  of  3  mm 
and  10%  active  metal  was  packed  into  the  bed.  Fig.  4(a-c) 
report  temperature  profiles  of  the  catalyst  and  the  outer  wall 
at  different  GHSVs.  From  the  results,  it  is  found  that  the 
temperature  difference  between  the  catalyst  and  outer  wall 
increases  markedly  as  GHSV  increases.  Heat  transfer  limi¬ 
tation  and  short  residence  time  degrade  the  performance  of 
the  SR  reactor  in  the  high  GHSV  region.  A  more  detailed 
discussion  of  this  behaviour  is  presented  below  in  Section 
3.6 
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3.  Mathematical  formulation 

3.1.  Governing  equations 

In  this  section,  mathematical  governing  equations  are  intro¬ 
duced  in  order  to  describe  the  physical  phenomena  in  the  steam 
reformer.  Mass,  momentum,  energy  and  species  equations  are 
solved  simultaneously  and  include  the  catalytic  reaction  on  the 
surface.  In  the  porous  medium,  permeability  is  a  key  parameter 
to  estimate  the  pressure  drop  in  the  flow.  For  a  more  rigorous 
analysis,  permeability  is  determined  by  the  lattice  Boltzmann 
method  [16]. 

For  the  chemical  reaction  on  the  surface  of  the  supported 
nickel  catalyst,  the  Langmuir-Hinshelwood  model  is  incorpo¬ 
rated,  assuming  that  this  reaction  is  the  rate-determining  step.  All 
the  governing  equations  are  formulated  in  an  axisymmetric  coor¬ 
dinates  system.  For  the  energy  equation,  a  heterogeneous  model 
of  gaseous  species  and  solid  catalyst  particles  is  employed  to 
consider  thermally  non-equilibrium  phenomena.  The  equations 
are  as  follows: 


dpf  3  13 

—  +  — (Pf^D)  +  -  —irpfV  D)  =  0 
at  dz  r  dr 

where  the  subscript  D  denotes  the  Darcian  velocity. 
dpfUB  3  .  1  3  (  . 

— X - h  —  (PfMDMD)  +  O'PfMDl’D) 

at  dz  r  dr 

dp  3  /  3wd\  13/  duv>\  fJ 

= - -|-  —  I  i  -  I  -I- - I  r  i  -  1  —  — i 

3  Z  dz 


(2) 


dz  J  r  dr 


dr  J  K 


(3a) 


3pf«D  ,3  %  ,  1  3 

— +  zr(PfUDVD)  +  -  —(mvDvD) 
dt  dz  r  dr 

dp  d  (  dvj)\  13/  dvj)\  pi 

=  +  fe  (WirJ  +  r*  (rMir)  -  Kn 


(3b) 


where  K  denotes  the  permeability. 

d(spfCv  f7f)  3  1  3 

aP’  +  T(PfCp,fuDTf)  +  -  —(rpfGpfVvTf) 
dt  dz  r  dr 

3  /  37f\  13/  dTf\  , 

—  T-  (  J  4  )  T  htat^Ts  —  Tf) 

dz  \  dz  J  r  dr  \  dr  J 


(4) 


where  &eff,f  denotes  the  effective  thermal  conductivity  for  bulk 
gas. 

3((1  ~  g)PsCp,sTs) 

dt 

3  /  3 Ts\  13/  3TS\  7 

—  (  ^eff,s^  )  H  zr  (  r^eff,s~  )  H-  htat,sf(Tf  —  Tf) 

dz  \  dz  J  r  dr  \  dr  J 


N 


+Pcat^( 

7=1 


(5) 


where  &eff,s  denotes  the  effective  thermal  conductivity  for  the 
catalyst. 

d(spfYf)  3  1  3 

+  -(puDYf)  +  --(rpvDYf) 

dt  dz  r  dr 

_  9  /n  9pfYf\  13/  3pfTf\ 

dz  V  e£f  dz  )  +  r  dr  V  Drff  dr  J 

T  Pf/?d^d,sf(^s  —  Yf) 

where  Y  is  the  mass  fraction  of  species. 

3((1  -  e)pfYs) 


(6) 


dt 


Pfhdad,sf(Yf  Fs)  +  pca.tri^i 


(7) 


where  r\  represents  the  conversion  rates  of  the  individual  species. 

3.2.  Boundary  conditions 

For  simplicity  of  the  given  problems,  the  transient  term  is 
omitted.  Hence,  only  the  boundary  conditions,  except  the  initial, 
condition,  are  implemented  to  obtain  the  solutions,  as  follows: 

(a)  At  the  reformer  inlet  face:  z  =  0 

37; 


7f  =  7f; 


dz 


=  0;  Yf  =  Yf 


(8) 


(b)  At  the  reformer  outlet  face:  z  =  L 

37f  dTs  3  If 

—  =0;  — =0;  — =0 

dz  dz  dz 

(c)  At  the  reformer  centre:  r  =  0 

—  =  0;  ^=0 

3r  3r 

(d)  At  the  interfacial  surface  of  the  inner  reformer  wall  and  the 


(9) 


(10) 


catalyst  bed:  r  =  R 


^eff,f  ^eff,S  —  <1  (z)»  7f  —  Ts 

dr  dr 

Tf=Ts  =  T(z ) 

3Ff 


dr 


=  0 


(Ha) 

(Hb) 

(He) 


Two  different  boundary  conditions  are  implemented  for  the 
energy  equation.  Eq.  ( 1 1  a)  is  used  for  the  heat  flux  condition,  and 
Eq.  (lib)  is  used  for  the  specified  reformer  wall  temperature. 

3.3.  Constitutive  equations 

In  order  to  solve  the  established  governing  equations,  some 
supplementary  relations  are  first  introduced  in  this  section.  Since 
the  SR  reaction  takes  place  in  porous  media,  all  the  physical 
parameters  should  be  modelled  in  this  region.  Permeability  can 
be  determined  by  the  following  correlation,  for  which  the  lattice 
Boltzmann  method  is  used,  as  follows: 

K  (  e11/3 

D2  V(l-e)2 


-c2] 


(12) 
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where  D  is  the  characteristic  length  of  packed  materials  and  the 
coefficients  C\  and  C2  are  dependent  upon  Kn0 ,  which  has  the 
following  correlations: 

Cl  =  0.709  -  l.62Kn0  +  5.892 Kn\  (13) 

C2  =  5.09  -  l4A4Kn0  +  36.84 Kn\  (14) 

where  Kn0  is  the  Knudsen  number  at  the  outlet. 

In  addition,  the  effective  thermal  conductivity  for  bulk  gas 
and  a  solid  can  be  determined  by 

*eff,f  =  *eff,s  =  (1  -  e)ks  (15) 

Also,  the  properties  of  gases  in  the  reformer  are  determined 
using  Chapman-Enskog  kinetic  theory  [17].  The  overall  prop¬ 
erties  such  as  density  and  heat  capacity  can  be  evaluated  based 
on  mole  or  mass  fractions  and  bulk  gas  temperature.  The  bulk 
gas  density  is  given  by 

Pf  =  - -  (16) 

RTElYi/Wi 

The  ideal  gas  law  is  incorporated  for  mixtures  of  gases. 

Since  a  heterogeneous  reaction  model  is  incorporated  in  the 
present  work,  the  heat  transfer  coefficient  and  interfacial  surface 
area  between  the  bulk  gas  and  catalyst  are  significant  deter¬ 
minants  of  reformer  performance.  The  heat  transfer  coefficient 
between  the  catalyst  and  the  gas  is  determined  by  the  Reynolds 
number  and  Prandtl  number.  Correlation  models  have  been  pre¬ 
sented  in  detail  in  the  literature  and  are  widely  accepted  [18]. 
The  interfacial  surface  area  between  the  bulk  gas  and  the  catalyst 
can  be  estimated  from  the  BET  surface  area  of  the  given  catalyst. 
Physical  variables  for  the  simulation  are  shown  in  Table  1 . 

3.4.  Chemical  reaction 

In  this  study,  SR  and  water-gas  shift  (WGS)  reaction  for 
methane  are  simultaneously  taken  into  account  to  analyze  the 
performance  of  the  reformer.  The  kinetic  rate  equations  for 
the  SR  and  WGS  reactions  are  taken  from  a  previous  study 
[19]. Three  overall  chemical  reactions  are  adopted: 

(i)  the  endothermic  SR  reaction  to  produce  H2  and  CO: 
CH4+H2O  **  CO  +  3H2 

Affi(298K)=  206,000kJkmol_1  (17) 

(ii)  the  WGS  reaction: 

CO  +  h2o  co2  +  h2 

AH1{ 298K)=  —41 ,  000 kJ kmol- 1  (18) 


Table  2 


Kinetic  parameters 


Reaction 

A{  (kmolkgca{  s  x) 

Ei  (kJ  mol-1) 

1 

1.174  x  1012bar°'5 

240.10 

2 

5.43  x  102  bar-1 

67.13 

3 

2.83  x  1011  bar0-5 

243.90 

(iii)  the  direct  SR  reaction  to  CO2: 

CH4  +  2H20  C02  +  3H2, 

A//3(298K)=  165,  000  kJ  kmol-1  (19) 


The  following  kinetic  rate  equations  are  derived  from  the 
Langmuir-Hinshelwood  model  assuming  the  surface  reaction 
is  the  slowest  step  among  adsorption,  surface  reaction  and  des¬ 
orption.  The  rate  equations  for  methane  [Eqs.  (17)— (19)]  can  be 
expressed  by 


n  =  L 

r2  =  k2 

=  k2 


Tch4Th2o/Ph2  ~  Pco^Hz/^pi 
DEN2 

PCOPHjO/ PH2  -  TC02/^p2 


DEN2 


PCH4/PH2  ~  PH1pc°JKPi 
DEN2 


DEN  is  defined  as 


(20) 

(21) 

(22) 


DEN  = 


1  +  KqoPco  +  Kn2pn2  +  ^ch4^ch4  +  Kh2oph2o 


Ph2 


(23) 


where  r,  (kmol  kg“J  h_1)  is  the  rate  of  reaction,  i.  Reaction  con¬ 
stants,  k\,  can  be  calculated  from  the  pre-exponential  factors 
using  the  Arrhenius  equation  and  the  Van’t  Hoff  equation.  The 
activation  energy,  Eu  and  adsorption  enthalpy  of  the  species, 
A77j,  are  taken  from  a  study  by  Xu  and  Froment  [19]: 


(24) 


K,  =  Aj  exp 


(25) 


Relevant  pre-exponential  factors,  activation  energy  and  equi¬ 
librium  constants  are  shown  in  Tables  2-4,  respectively. 


Table  1 

Physical  parameters  for  the  simulation 

Table  3 

Adsorption  constants 

Physical  variables 

Value 

Species 

Ai  (bar  1) 

AHi  (kJ  mol-1) 

Inlet  pressure  (pm )  (bar) 

1 

ch4 

6.65  x  nr4 

-38.28 

Porosity  (s) 

0.6 

CO 

8.23  x  lO"5 

-70.65 

Interfacial  surface  area  (as f)  (1  m-1 ) 

2000 

H2 

6.12  x  nr9 

-82.90 

Catalyst  density  (pCat)  (kgm-3) 

2000 

h2o 

1.77  x  lO"5  bar 

88.68 
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Table  4 

Equilibrium  constants 


Reaction 

Equilibrium  constant  7Tp; 

Unit 

1 

5.75  x  1012exp(— 11476/7) 

bar2 

2 

1.26  x  10-2  exp(4639/7) 

bar0 

3 

KpiKp2 

bar2 

3.5.  Code  validation 
3.5.1.  Cold  flow 

In  this  section,  a  developed  code  is  validated  with  analytical 
solutions  for  micro-channel  flow  as  a  benchmark  [20].  The  SR 
reaction  takes  place  in  a  packed  bed,  which  is  modelled  as  porous 
media.  Analysis  of  fluid  flow  in  porous  media  must  be  performed 
carefully,  since  the  pressure  drop  is  comparably  greater  than  in 
the  regular  fluid  region.  A  high  pressure  drop  affects  the  system 
efficiency  of  the  steam  reformer.  As  shown  in  Fig.  5(a),  the  ideal 
composite  system  consists  of  a  porous  medium  and  an  overly¬ 
ing  fluid  layer.  When  the  flow  enters  the  ideal  composite  system, 
the  normalized  velocity  profile  is  compared  with  the  analytical 
solution  in  the  developed  flow  region,  as  presented  in  Fig.  5(b). 
The  Darcy  number  in  the  simulation  is  0.001,  and  continuous 


(b) 


(b)  Furnace  temperature(°C) 


Fig.  6.  (a)  Comparison  of  temperature  profiles  between  numerical  and  experi¬ 
mental  results  vs.  different  GHSVs  (7™  =  480  °C,  S:C  =  3.0).  (b)  Comparison  of 
selected  gas  compositions  between  numerical,  thermodynamic,  and  experimen¬ 
tal  results  vs.  different  furnace  temperatures  (7™  =480  °C,  S:C  =  3.0). 


velocity  and  slip  shear  stress  conditions  are  imposed  at  the  inter¬ 
face  between  the  regular  fluid  and  the  porous  media.  From  the 
results,  it  can  be  verified  that  the  developed  code  predicts  the 
fluid  flow  in  the  porous  media  with  adequate  accuracy,  as  shown 
in  Fig.  5(b). 

3.5.2.  Reacting  flow 

Using  the  established  model,  numerical  simulations  are  car¬ 
ried  out  (using  the  same  geometry  and  operating  conditions 
as  applied  in  the  experiment)  to  validate  the  developed  code. 
Temperature  profiles  measured  from  the  outer  reformer  wall  are 
implemented  as  a  thermal  boundary  condition  to  compare  the 
results.  The  data  given  in  Fig.  6(a)  and  (b)  conform  that  the 
numerical  results  are  in  good  agreement  with  the  experimental 
results. 

3.6.  Numerical  results 


Fig.  5.  (a)  Ideal  composite  system;  and  (b)  comparison  between  numerical  and  this  section,  numerical  investigations  are  made  of  the 

analytical  solution  ( s  =  0.7).  effects  of  parameters  for  operating  conditions  such  as  inlet  tern- 
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Fig.  7.  Bulk  gas  temperature  and  methane  conversion  vs.  different  inlet  gas 
temperatures  (7W  =  800°C,  S:C  =  3.0,  GHSV  =  1000 h-1). 


Fig.  8.  Bulk  gas  temperature  and  methane  conversion  vs.  different  reformer  wall 
temperatures  (7in  =  450°C,  S:C  =  3.0,  GHSV  =  3000  h”1). 


perature,  reformer  wall  temperature,  S:C  ratio,  and  GHSV.  It 
is  found  that  high  values  of  inlet  gas  temperature  and  reformer 
wall  temperature  are  favourable.  Nevertheless,  it  is  necessary  to 
take  energy  efficiency  into  account  at  high  operating  tempera¬ 
tures.  To  investigate  the  SR  reformer  numerically,  a  bench-scale 
catalyst  bed  of  6  cm  in  diameter  and  30  cm  in  length  was 
selected. 

To  assess  the  influence  of  the  inlet  gas  temperature,  the  reactor 
wall  is  maintained  at  800  °C  the  S:C  ratio  at  3.0,  and  the  GHSV 
at  1000  h-1.  When  the  inlet  gas  temperature  is  higher,  the  cat¬ 
alyst  bed  temperature  at  the  centre  is  also  higher  (see  Fig.  7), 
which  results  in  higher  production  of  hydrogen  and  methane 
conversion.  Methane  conversion  is  defined  as 


^CH4  -  *CK 4(Z*) 


Z7  in 
^CH4 


(26) 


where  is  the  flux  of  methane  at  the  inlet,  and  Fch4(z*)  is 
the  flux  of  methane  at  z  =  z*. 

In  Fig.  7,  the  catalyst  bed  can  be  classified  into  three 
regions — (1)  region  1:  reaction-dominant  region;  (2)  region  2: 
mixed  reaction  and  heat  transfer  region;  and  (3)  region  3:  heat- 
transfer-dominant  region.  At  the  entrance  region  (region  1),  the 
gas  mixture  enters  and  SR  predominantly  takes  place,  resulting 
in  a  decrease  in  the  bulk  gas  temperature.  The  heat  transfer  and 
SR  then  occur  simultaneously  in  region  2.  Here,  the  heat  trans¬ 
fer  is  marginally  greater  than  the  reaction  and  therefore  the  bulk 
gas  temperature  continuously  increases.  Afterwards,  in  region 
3,  heat  transfer  from  the  wall  is  dominant,  since  there  is  not 
sufficient  fuel  to  reform,  so  that  the  gradient  of  the  temperature 
profile  is  steeper  than  that  of  region  2. 

The  influence  of  the  reactor  wall  temperature  was  also  inves¬ 
tigated.  The  inlet  gas  temperature  is  450  °C,  the  S:C  ratio  is  3.0, 
and  GHSV  is  3000  h_1 ,  which  constitutes  a  relatively  short  res¬ 
idence  time  with  this  geometry.  As  an  operating  parameter,  a 
maximally  high  reactor  wall  temperature  is  also  favourable  with 
respect  to  energy  efficiency.  The  data  in  Fig.  8  show  that  the  bulk 
gas  temperature  behaviour  can  be  clearly  distinguished  on  the 


Relative  lengh  (z/L) 


Fig.  9.  Bulk  gas  temperature  and  methane  conversion  vs.  different  S:C 
(r1  =450°C,  Tw  =  800°C,  GHSV=  lOOOh-1). 


Fig.  10.  Bulk  gas  temperature  and  methane  conversion  vs.  different  GHSVs 
(Tn  =  450  °C,  Fw  =  800  °C,  S:C  =  3.0). 


Methane  conversion  §  Methane  conversion 


Methane  conversion  Bulk  gas  temperature, Tc(°C) 
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(b) 


Fig.  11.  (a)  Typical  packing;  and  (b)  mixed  packing. 


basis  of  two  regions;  region  3  as  shown  in  Fig.  7  does  not  arise 
due  to  the  relatively  higher  GHSV. 

The  effect  of  S:C  ratio  when  the  inlet  gas  temperature  and 
reformer  wall  temperature  are  maintained  at  450  and  800  °C, 
respectively,  is  presented  in  Fig.  9.  The  GHSV  is  1000  h-1.  As 
the  S:C  ratio  is  increased,  methane  conversion  also  increases. 
An  overly  high  S:C  ratio  can,  however,  lead  to  relatively 
low  hydrogen  concentrations  according  to  the  chemical  equi¬ 
librium.  Hence,  an  appropriate  range  of  S:C  ratio  must  be 
selected. 

Finally,  the  influence  of  GHSV  was  examined  based  on 
the  heat  transfer  limitation  from  the  reactor  wall.  The  inlet 
gas  temperature  and  wall  temperature  were  450  and  800  °C, 
respectively,  and  the  S:C  ratio  was  3.0.  Methane  conversion  is 
found  to  be  remarkably  poor  when  GHSV  exceeds  5000  h-1, 
see  Fig.  10.  Hence,  the  GHSV  should  be  carefully  selected  as 
regards  the  residence  time  of  the  gas  mixture  and  the  effect 
of  heat  transfer  from  the  wall  to  the  centre  of  a  given  catalyst 
bed. 


4.  New  method  of  catalyst  packing 


0.0  0.2  0.4  0.6  0.8  1.0 

Relative  length  (z/L) 

Fig.  12.  (a)  Bulk  gas  temperature  vs.  different  packing  methods  and  GHSVs 
(heat  flux:  5kWm-2,  7in  =  450°C,  S:C  =  3.0);  and  (b)  methane  conversion 
vs.  different  packing  methods  and  GHSVs  (heat  flux:  5  kW m-2,  7™  =450  °C, 
S:C  =  3.0). 


In  this  section,  an  efficient  catalyst  packing  method  is  pro¬ 
posed  with  the  aim  of  enhancing  the  heat  transfer  characteristics. 


Fig.  13.  (a)  Bulk  gas  temperature  vs.  different  packing  methods  and  GHSVs 
(heat  flux:  lOkWm-2,  7in  =  450°C,  S:C  =  3.0);  and  (b)  methane  conversion 
vs.  different  packing  methods  and  GHSVs  (heat  flux:  10 kW  m-2,  Pn  =  450  °C, 
S:C  =  3.0). 
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Fig.  14.  (a)  Bulk  gas  temperature  vs.  different  packing  methods  and  GHSVs 
(heat  flux:  15kWm-2,  7in  =  450°C,  S:C  =  3.0);  and  (b)  methane  conversion 
vs.  different  packing  methods  and  GHSVs  (heat  flux:  15  kWm-2,  Pn  =450  °C, 
S:C  =  3.0). 


Fig.  15.  Hydrogen  production  to  fed  methane  (mole  to  mole)  vs.  different  pack¬ 
ing  methods,  GHSVs,  and  heat  fluxes  (Tn  =  450  °C,  S:C  =  3.0). 


Fig.  16.  Maximum  inner  tube  wall  temperature  vs.  different  packing  methods, 
GHSVs,  and  heat  fluxes  (Fn  =  450  °C,  S:C  =  3.0). 


Two  different  catalyst  packing  configurations  are  illustrated  in 
Fig.  11,  namely,  (a)  a  conventional  packing  and  (b)  new  config¬ 
uration  to  enhance  heat  transfer  and  the  endothermic  reaction. 

In  the  new  packing  configuration,  inert  catalysts  are  inserted 
between  active  catalysts  in  a  series  to  recover  heat  at  the  region 
where  heat  transfer  from  the  wall  occurs  and  there  is  no  signif¬ 
icant  endothermic  reaction.  Thus  the  gas  temperature  increases 
at  this  region  and  potentially  enhances  the  desired  reaction  at  the 
next  catalytic  region.  The  performance  of  this  packing  method 
has  been  extensively  compared  with  typical  methods  under  the 
same  operating  conditions.  In  the  new  packing  approach,  cata¬ 
lyst  loading  is  reduced  by  50%,  which  reduces  the  production 
cost  of  the  syngas. 

In  Figs.  12-14,  (a)  illustrates  the  bulk  gas  temperature  at 
the  centre  of  the  reactor,  whereas  (b)  gives  methane  conversion 
with  different  heat  fluxes  and  GHSVs.  In  Fig.  12(a)  and  (b),  the 
temperature  at  the  reactor  centre  and  the  methane  conversion 
are  compared  when  a  heat  flux  of  5  kW  m-2  is  applied.  Methane 
conversion  and  hydrogen  production  are  similar  at  both  GHSVs, 
as  shown  in  Figs.  12(a)  and  15.  Nevertheless,  the  temperature  of 
the  inner  tube  wall  with  the  new  packing  approach  is  maintained 
at  a  slightly  higher  level  than  that  experienced  with  a  typical 
packing  approach,  as  shown  in  Fig.  16. 

In  Fig.  13(a)  and  (b),  more  heat  is  supplied  to  the  reactor.  The 
hydrogen  production  of  the  proposed  packing  configuration  is 
similar  to  that  of  a  typical  packing  approach,  as  shown  in  Fig.  15. 
The  maximum  inner  tube  wall  temperature  with  the  new  pack¬ 
ing  method  is  remarkably  higher  at  the  rear  of  the  bed,  because 
the  endothermic  reaction  is  finished  at  a  relatively  low  GHSV 
(2000  h-1).  In  this  case,  a  slightly  higher  GHSV  (5000  h-1) 
reduces  the  tube  wall  temperature  (note,  the  GHSV  should  be 
also  selected  with  care  in  view  of  hydrogen  production).  Hydro¬ 
gen  production  with  respect  to  the  supplied  methane  at  the  same 
GHSV  is  similar  to  that  for  a  typical  packing  configuration,  as 
shown  in  Fig.  15. 

The  bulk  gas  temperature  and  methane  conversion  when  heat 
flux  of  15  kW  m-2  is  applied  are  shown  in  Fig.  14(a)  and  (b), 
respectively.  For  a  relatively  low  GHSV  (2000  h-1,  case  3), 
hydrogen  production  with  the  new  packing  method  is  better  than 
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that  obtained  with  the  typical  packing  method.  In  the  latter,  the 
heat  activates  the  endothermic  reaction  very  strongly  near  the 
reactor  wall,  compared  with  the  other  cases  (cases  1  and  2). 
Thus  heat  cannot  penetrate  to  the  centre  of  the  reactor  and  this 
results  in  low  methane  conversion  and  hydrogen  production.  On 
the  other  hand,  in  the  relatively  high  GHSV  (5000  h-1)  case, 
hydrogen  production  is  almost  the  same  for  the  different  pack¬ 
ing  methods.  These  phenomena  result  from  the  relatively  low 
residence  time  for  the  reactions  in  the  reformer.  Nevertheless, 
the  mixed  packing  method  is  still  advantageous  for  economical 
use  of  catalyst. 

Maximum  inner  tube  wall  temperatures  for  the  two  different 
packing  methods  are  given  in  Fig.  16.  When  a  relatively  small 
amount  of  heat  is  applied  (case  1),  the  maximum  temperature 
difference  between  the  two  packing  methods  is  quite  small  so 
that  thermal  impact  on  the  reformer  tube  is  virtually  the  same. 
On  the  other  hand,  material  problems  could  arise  in  the  high 
heat  flux  region  (cases  2  and  3),  as  shown  in  Fig.  16.  These 
problems  can  be  solved  by  the  operating  in  a  relatively  higher 
GHSV  region  (i.e.,  5000  h-1  compared  with  2000  h-1). 

5.  Conclusion 

Steam  reforming  of  methane  over  nickel  alumina  has  been 
investigated  both  experimentally  and  numerically.  Based  on  the 
experimental  results,  a  heterogeneous  two-dimensional  reac¬ 
tor  model  has  been  developed  and  extensively  validated.  The 
influence  of  the  operating  parameters  (e.g.,  inlet  gas  temper¬ 
ature,  reactor  wall  temperature  S:C  ratio  and  the  GHSV)  has 
been  examined.  It  is  found  that  high  inlet  gas  and  reactor  wall 
temperatures  are  favourable,  but  these  parameters  have  to  be 
selected  carefully  with  respect  to  energy  efficiency.  An  appro¬ 
priate  S:C  ratio  should  be  selected  in  order  to  maximize  the 
hydrogen  production  rate.  On  the  other  hand,  a  high  gas  flow 
rate  is  unfavourable  due  to  the  short  residence  time  and  the 
heat  transfer  limitation.  Thus,  a  moderate  GHSV  also  should 
be  applied.  From  these  results,  it  is  concluded  that  heat  trans¬ 
fer  is  very  significant  to  obtain  higher  hydrogen  production  in  a 
given  catalyst  bed. 

To  overcome  the  heat  transfer  limitation  at  moderate  and 
higher  flow  rates,  and  thereby  improve  reactor  performance,  a 
new  catalyst  packing  method  is  proposed.  When  heat  is  applied 
to  the  reactor,  the  hydrogen  production  is  almost  the  same  or 
better  than  that  obtained  with  a  typical  packing  approach.  The 


proposed  catalyst  packing  method  shows  good  performance  and 
requires  a  smaller  loading  of  catalyst  because  heat  recovery  is 
carried  out  in  an  inert  catalyst  bed.  It  should  be  noted,  however, 
that  the  reactor  tube  for  the  mixed  packing  method  should  be 
carefully  designed  so  as  to  avoid  hot  spots. 
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